#include "evaluar.h"
#include "spline.h"

int main(int argc, char** argv)
{
	int n = 6;
	
	vPunto puntos;
	
	forn(i, n)
	{
		puntos.push_back(Punto(i*i, i*i+2*i));		
	}
	
	forn(i, n)
		cout << "(" << puntos[i].x << ", " << puntos[i].y << ")" << endl;
	
	vSpline spline (n-1);
	
	interpolar(puntos, spline, n);
	
	forn(i, n-1)
	{		
		double punto = spline[i].x_i;
		cout << "S_" << i << "(x) = " << spline[i].d << "(x-" << punto << ")^3 + ";
		cout << spline[i].c << "(x-" << punto << ")^2 + ";
		cout << spline[i].b << "(x-" << punto << ") + ";
		cout << spline[i].a << endl;
	}
	
	cout << "Verificando spline..." << verificar((vPunto&) puntos, (vSpline&) spline, n) << endl;
	
	return 0;
}
